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ABSTRACT 



Of the fourteen transiting extrasolar planetary systems for which radii have been 
measured, at least three appear to be considerably larger than theoretical estimates 
suggest. It has been proposed by Bodenheimer, Lin & Mardling that undetected com- 
panions acting to excite the orbital eccentricity are responsible for these oversized 
planets, as they find new equilibrium radii in response to being tidally heated. In 
the case of HD 209458, this hypothesis has been rejected by some authors because 
there is no sign of such a companion at the 5 ms""'^ level, and because it is difficult 
to say conclusively that the eccentricity is non-zero. Transit timing analysis as well 
as a direct transit search have further constrained the existence of very short-period 
companions, especially in resonant orbits. Whether or not a companion is responsible 
for the large radius of HD 209458b, almost certainly some short-period systems have 
companions which force their eccentricities to nonzero values. This paper is dedicated 
to quantifying this effect. 

The eccentricity of a short-period planet will only be excited as long as its (non- 
resonant) companion's eccentricity is non-zero. Here we show that the latter decays on 
a timescale which depends on the structure of the interior planet, a timescale which 
is often shorter than the lifetime of the system. This includes Earth-mass planets in 
the habitable zones of some stars. We determine which configurations are capable of 
sustaining significant eccentricity for at least the age of the system, and show that 
these include systems with companion masses as low as a fraction of an Earth mass. 
The orbital parameters of such companions are consistent with recent calculations 
which show that the migration process can induce the formation of low mass planets 
external to the orbits of hot Jupiters. Systems with inflated planets are therefore good 
targets in the search for terrestrial planets. 

Key words: planetary systems - celestial mechanics - stellar dynamics - methods: 
analytical - planetary systems: formation 



1 INTRODUCTION 



While most extrasolar planet discoveries to date have relied on the radial velocity method of detection with its inherent 
uncertainty in the planetary mass, a few systems have been observed in transit, thereby breaking this degeneracy. Not only 
does this allow a reasonably accurate determination of the planetary mass (given a reliable estimate for the stellar mass is 
available), it allows on e to estimate the radius of the planet, and hence its density. Before the discovery that HD 209458b 
transited its host star ( Charbonneau et al. 200d : Henrv et al. 2000l l . it was by no means clear that the planets discovered till 
then, all with masses similar to or more than Jupiter's, were gas giants. The extra information gleaned from photometric transit 
data has provided the opportunity to stud y planetary structure JSodenheimer. Lin fc MardlinglboOl ). while spectroscopic data 
have revealed their chemical composition IjCharbonneau et al.ll2002l ). In the former case, determination of the planetary radius 
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Table 1. Observational data for inflated planets 



Planet 


Period (day) 


(Mq) 


nip {Mj) 


ap (AU) 




Rp (Rj) 


HD 209458b 
WASP-lb_^ 
HAT-P-lb 


3.5247 
2.5199 
4.4653 


1.10 
1.15 
1.12 


0.64 ± 0.06 

79+0.13 
'-'•'^-0.06 

0.53 ± 0.04 


0.045 
0.0379 
0.0551 


0.014 ± 0.009^ 
0.09 ± 0.02_^ 


1.32 ±0.03 

1.443 ± 0.039 

1 36+°-" 
^•'-'"-0.09 



iLaughlin et~al] l|2005a). 
° IStempels et al. 

^ This value is based on very little data l lBakos et al.ll2006l) . 



allows one to put constraints on the possible existence of a planetary core, and this in turn allows one to constrain the mode 
of planet formation, at least for that planet. In pri nciple this allows a choice between the c ore accretion model |Sa fronov 



196C 



Wetherii] 1980l : Bodenheimer fc Pollackl 19861 ) and the gravitational instability model ( Kuiper 195ll : CameronF 



1978 



Bosslll997l ). although it is not clear that the latter should produce core-less planets given accretion of solids plus subsequent 



gravitational settling almost certainly follows formation. Either way, such an exercise requires knowledge of the equilibrium 
structure of a gas giant planet, with and without a core, and with and without external and internal heating. 

While it can be argued that there are problems matching theoretical radii with observation for all transiting planets 
j Burrows et al.ll2007 ) , thre e stand out a s particularly excessive. Thes e are HD 209458b ( Charbonneau et allboool : Henry et al 



200' 



Knutson et al 



20o3), WASP- lb (jCoUier Cameron et ahlbooel l and HAT-P-lb jBakos et al.lbood '). for which the plan 



etary radii and other observationally derived parameters are listed in Table [T] Data for HD 209458 and HAT-P-1 are taken 
from Burrows et al. (2006) and references therein, except for the eccentricities, while data for WASP-1 are taken from 
Stempels et al.l l|2007l ). Here m, and nip are the masses of the star and planet respectively, ap and Bp are the orbital semima- 
jor a xis and eccentri c ity and Rp is the estimated planetary radius {Rj is Jupiter's radius). The eccentricity of HD 209458b is 
from Laughlin et al. |2005al). a paper devoted specifically to discussing this contentious orbital element, while that of HAT-P- 



200a), a value based on very little data. However, as pointed out by the latter authors, HAT-P-lb is 
further from its parent star than HD 209458b, and is such that the relatively large estimated value of 0.09 produces a similar 
planetary radius. 

Various mechanis ms have been suggested to account for the excess radii of the planets listed in Table U Obvious contend ers 
are stellar irradiation ( Burrows et al ■hoodlGuillot fc Showmaril2002l:lBodenheimer Laughhn fc Linl200a :, Burrows. Sudarskv fc HubbardI 



20031: lArras fc Bildstenlbood: iFortney Marlev fc Barnesll2007l: I Burrows et al.l 120071 ) and tidal forcing by a companion planet 



( Bodenheimer. Lin fc Mardlin j " 2001 ; Bodenheimer. Laughlin fc Lin 20031 ). the latter requiring non-zero eccentricity of the 



observed planetary orbit to be viable. 

A less obvious mechanism has been proposed by IWinn fc Holman (I2OO5I) in which "o bliquity tides" are forced if the 
system can manage to remain in a suitable Cassini state. Both iLevrard et al.l (12007|) and iFabrvckv. Johnson fc Goodman 
1 20071 ) verify that in principle it is possible for this mechanism to operate as long as the system dissipates energy at a rate 
which is too small to inflate the planet to the observed size. 

This paper focuses on the mechanism of tidal forcing provided by a companion planet, and hence requires that the orbit 
of the transiting planet be non-circular. Us i ng a bo otstrap Monte Carlo technique to estimate confidence intervals for the 



eccentricity of HD 209458b, iLaughhn et al.1 (|2005al ) find that Bp < 0.049 with 99.73% confidence, Bp < 0.034 with 95.4% 



confidence, and Cp < 0.02 with 68.3% confidence, with longitude of periastron ro^ = 80° ± 80°. Their best-fit measurement 
for the eccentricity is Cp — 0.014. However, they point out that because ep has a hard lower bound of zero, any best-fit model 
to a radial velocity data set will have ^ 0, with largest estimates corresponding to vop ~ 90° and 270°. Using synthetic 
data corresponding to a truly circular orbit, such bootstrap trials produce a distribution in vop — ep space which is symmetric 
about these two values, while use of the real radial velocity data produces a distinctly asym metric distribution whic h is biased 
towards vop = 90°. This suggests that in fact 7^ 0, although gp = cannot be ruled out. Laughlin et al. ( 2005a[ ) point out 



that the fact that the estimate for zop is centred near 90° suggests that even if HD 209458b has a no n-zero eccentricity, timing 
of the secondary eclipse will find cos Wp ~ 



as has recently been found bv lDeming et al.l (j2005l ) in the infrared. Note also 
that the orbital model of ILaughlin et al.l ( 2005al ) leads to a stellar jitter estimate of 3.3 ms~^ which is close to the predicted 
value of 2.8 ms"^ Thus, at least for HD 209458, we concentrate our attention in this paper on companions which produce 
stellar reflex velocities not much larger than this value. 

The presence of a very short-period companion to HD 209458b has been ruled out by transit-ti ming analysis for inte rior 
and exterior planets with periods up to 21 days (corresponding to the 6:1 resonance at 0.16 AU ) jAgo l fc Steffen**2007l). as 
well as by a direct transit search for planets with radii down to 2_R® and periods up to 14 days (|Croll et al.ii2007l ). Here we 
will concentrate on companion planets with semimajor axes greater than these lower bounds. 

This paper presents a generic analysis for short-period systems with external companion planets (low-mass or otherwise). 
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estimating the range of orbital parameters and masses required for such companions to be responsible for inflated radii. An 
unexpected result is that very low mass planets are sometimes capable of inducing significant tidal forcing on extremely 
long timescales. Section 2 discusses the relationship between t he eccentricity and the Q -value of a planet for a given value 
of the planet's luminosity, the Q- value or tidal quality factor ( Goldreich fc Soter 19661 ) being a measure of how efficiently 
the tides dissipate their o scilla tion energy. Section 3 discu sses three-body dynamics in the presence of dissipation (see also 
Murray fc Dermotj (|2000l ') and I Adams fc Lauehli nl (|2006bh ). We start by review ng the secular evolution of point-mass copla- 



nar three-body systems, with and without relativistic effects, deriving explicit formulae for the periods of variation of the 
eccentricities and their amplitudes. We identify the orbital parameters which correspond to libration and circulation of the 
angle 77 = u^p — zxjc, where -ooc is the longitude of periastron of the companion, and in particular find the fixed points in ep — 77 
space. We then show that tidal damping results in the system evolving towards the relevant fixed point (now a quasi-fixed 
point) on three times the tidal damping timescale. Once this is achieved the system continues to evolve toward the true fixed 
point for which both eccentricities are zero. The timescale for this latter phase may be several orders of magnitude longer 
than the tidal damping timescale, and can in fact be much longer than the age of the system. We illustrate the theory for 
hypothetical companions to HD 209458, and show that companion masses as low as a fraction of an Earth mass are capable 
of sustaining non-zero eccentricites in the observed planetary orbit for at least the age of the system, and in some cases at 
the level required to sustain the observed radii at the current epoch. 

Section 4 maps out regions in rric — ac parameter space, where rric and are the mass and semimajor axis of a hypothetical 
companion, for which the eccentricities of HD 209458b and HAT-Plb would be at their observed level, as well as for WASP-lb 
which for which the eccentricity is currently unconstrained. Section 5 applies the secular theory developed here to determine 
which configurations have circularized companion planets, while Section 6 presents a summary and a discussion. 



2 THE ECCENTRICITY-RADIUS RELATION 



Various authors have studied the problem of the response of a planet to internal and external heat sources (eg. Burrows et 
al. 2000, Bodenheimer, Lin fc Mardling 2001, Bodenheimer, Laughlin fc Lin 2003, Burrows et al. 2007). In particular, models 
include planets with and without cores, with and without stellar irradiation, and with and without an internal heat source. In 
all these models the planet assumes an equilibrium radius, Rp, and it is the relationship between this and the rate of energy 
dissipation in the planet, Cp, that is used in an attempt to explain the range of planetary radii observed in transiting planets. 

In the case that jCp is due to tide s raised by an external companion, the rate of internal energy dissipation in the planet 
IS given by feg.. lYoder fc Pealel lfl98lh . lBodenheimer. Laughlin fc LinI (|2003h ) 



Gmtrrir. 



p 

Tcirc 



to second order in the eccentricity, where 

5 



21ni, 



Op \ 

RpJ 



(1) 



(2) 



is the circularization timescale of the orbit, with Up the orbital frequency, kp the tidal Love number of the planet and Qp its 
Q-value or tidal quality factor (Goldreich fc Soter 1966)0 Equating this to the power, £p, needed to maintain an observed 
planetary radius gives Bp in terms of the observable quantities m*, Rp and Op {mp actually drops out), as well as the theoretical 
quantities Qp, kp and £p: 



Mp, [RjJ V0.04AUy \10^J \kj 



-1/2 



10-8L. 



1/2 



(3) 



where kj = 0.34 is the tidal Love number for Jupiter. iBurrows et al.l (|2007l ) calculate £p for c areless solar metalicity planets 
which are irradiated by their parent stars, values for which are listed in Table [2] under Cp^ . I Boden heimer. Laughlin fc Lin 



( 20031 ) calculate £,p specifically for HD 209458b only (in addition to various generic models), and find a value for Cp ten times 
higher for a coreless planet, while for a planet with a 20 Af® core, they find £p is another factor of ten higher. Note that they 
take a slightly higher value for Rp {1.35 Rj). Values for Cp corresponding to Cp — Cp^ and Qp = 10^ (see discussion below) 
are shown in Table [2] as well as values corresponding to Cp = \Q Cp^ . The estimated value of Cp for HD 209458b (listed under 
e°l"') is more consistent with Cp — Cp^ than Cp = 10 Cp" while the opposite is true for HAT-P-lb, although the eccentricity 



^ Note that Goldreich & So ter (1966) use a modified Q-value, Qp, which absorbs the Love number such that Qp = 3Qp/2fcp; equation 
(4) in (Bodenheimer, Laug hUn &: Linl l200^ should have Qp instead of Qp. Another potential source of confusion is the use of the 
(quadrupole) apsidal motion constant k2 rather than the tidal Love number; these are related by kp = 2k2. 
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Table 2. Theoretical data assuming Qp = 10^ 



Planet 


Ago (Gyr) 


rcirc (Gyr) 


{Lq) 




ep{10£S") 


^ohs 


HD 209458b 


5.5 ±1.5^ 


0.045 


1.02 


X 10"* 


0.014 


0.046 


0.014 ±0.009 


WASP- lb 


2.0 ± 1.0^ 


0.011 


4.49 


X 10"* 


0.012 


0.037 




HAT-P-lb 


3.6±1.0~ 


0.109 


1.32 


X 10-8 


0.031 


0.098 


0.09 ±0.02 



" iBurrows et al] l[2003) 
^ IStempels et al. I l l2007^ 



value for the latter awaits refinement. Note that in reality Cp is a function of Rp and nip, however, its role in equation ^ is 
to parameterize various a uthors' models. 

In both the models of Burrows et al. (2007) and Bodenheimer. Laughlin &: Lin ( 20031 ) . tidal energy is deposited homoge- 
neously throughout the planet's envelope and dissipated locally. If, for example, tidal energy is predominantly dissipated in 
the outer layers of the planet, a significantly smaller value of Cp, and hence e,, , may be r equired to achieve the same planetary 
radius (as long as the energy is dissipated sufficiently deep below the surface; IWu finds that tidal excitation of inertial 

modes results in energy being dissipated very close to the surface resulting in very little structural adjustment). Alternatively, 
for a given ep one would require larger values of Qp for correspondingly smaller values of Cp. Note that this scenario makes 
it difficult to draw conclusions about the existence or otherwise of a planetary core. 

There is significant uncertainty in the theory of tidal damping in stars and planet s, much of it associa ted with internal 
structure and (in the case of gaseous objects) the theory of turbulent dissipation (see Qgilvie fc Lin ( 2004h for an excellent 
review). In the case of solid planets, one can infer damping timescales based on laboratory measurements of Q-values of solid 
materials, and for the special case of the Earth-Moon system, one can actually measure the rate of recession of the Moon. In 
the case of Jupiter and its satellite system, it is possible to put lower and upper bounds on Jupiter's Q-value (6 x 10* - 2 x 10*^ : 
Yoder & Peale 1981). These bounds are inferred from the existence of the Laplace resonance: if Qj were too low lo would 
have already passed through the 2:1 resonance with Europa (probably destabilizing th e system), while if it were t oo high it 
would never have been captured in the first place. Similarly for Uranus, Qu < 39, 000 jTittemore fc Wisdomlll989l ) while for 
Neptune, 12,000 <Qn < 330,000 (jBanfield fc Murravl[l99a ). 

There are various factors which might be expected to influence the Q-value of a hot Jupiter and how this might differ 
from that of Jupiter. One obvious difference between hot Jupiters and our own Jupiter are spin rates. While Jupiter's spin 
period is 10 hrs, hot Jupiters are approximately synchronous with the orbital motion and hence have periods around an order 
of magnitude longer. Goldreich fc Soter (1 966, p388) poin t out that the smaller the mass of a planet, the lower the Q-value. 
More recent work by Qgilvie fc Lin ( 20041 ) and Wul ( 2005 ) who consider tidal forcing in slowly rotating planets (approp riate 



to syn chronously rotating short-period planets) find that Q-values for such planets tend to be similar to that of Jupiter. IWu 



1 20031 ) shows that the observed orbital parameters of short-period extrasolar planets indicate that their Q-values lie between 
3 X 10* and 3 x 10^. The lower value is increased if one takes into account the fact t hat planets are likely to arrive at the ir 
current positions with significantly larger radii that those observed today (see Fig. 1 of Bodenheimer. Lin fc Mardlinj ( 2001 )). 
Throughout this paper we use mostly Qp — 10^ to illustrate the theory developed here. 



3 LONG-TERM ECCENTRICITY EVOLUTION 

The estimated ages of HD 209459, HAT-P-1 and WASP-1 are shown in Table [2] together with the circularization timescales 
of their planets (assuming Qp = 10^). Clearly a mechanism is needed to maintain orbital eccentricity if tidal heating is the 
source of energy which maintains the inflated planetar y radii. It seems to be imp licitly assumed that the oscillatory secular 



rnpl: 

variations in the eccentricities of a pair of close planets ijMurrav fc Dermottll200d ) prevail indeflnitely, even in the presence of 



tidal dissipation. In this section we show that the natural endpoint of such evolution is in fact a pair of circular orbits, with 
the structural properties of the tidally forced planet determining (at least in part) the timescale on which this occurs. This 
evolution occurs in three stages (see Figure [3]): (1) circulation of the angle between the lines of apsides of the two planets 
accompanied by the slow oscillation of the eccentricities at constant amplitude, together with decline of the mean value of 
the inner planet's eccentricity until it reaches a (quasi)-fixed value: this occurs on the circularization timescale; (2) libration 
of the angle between the lines of apsides accompanied by the slow oscillation of the eccentricities with reducing amplitude 
but maintenance of the mean value of the inner eccentricity: this occurs on twice the circularization timescale; and (3) either 
(a): a slow non-oscillatory decline in both eccentricities to zero or (b): maintenance of the inner eccentricity for some (often 
extremely long) period of time followed by its increase to some maximum value followed by an often fast decline of both 
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eccentricities to zero. Whether or not the final phase occurs via route (a) or (b), it proceeds on a timescale sometimes several 
orders of magnitude longer than the circularization timescale. 

In general, the rates of change of a planet's orbital elements will have contributions from companions if they exist, 
the relativistic post-Newtonian contribution to the potential of the star, the tidal and spin bulges of the star and planet, 



and tidal dissipation in the atmospheres of the star and the planet ( Mardling fc LinI 2002 ; Adams fc Lauehlinl 2006a 3) . For 



coplanar systems, planetary elements which evolve are the semimajor axis, the eccentricity and the longitude of periastron. 
The only effect which contributes to a secular change in the semimajor ax;is (except when the system is in a resonance) 
is tidal dissipation, while for short-period systems the apsidal motion is dominated by the third body and the relativistic 
potential of the star. The eccentricity in turn is dominated by the third body and tidal damping. We will assume that the 
orbits are coplanar, that the planet is synchronously rotating with the orbital motion, and that the spin bulge of the star is 
negligible. For the analysis we will assu me that the rate of chan ge of the planet's semimajor axis is negligible, while numerical 
integrations will include this effect (see Rasio et al. ( 1996h and Mardling fc Lin (2004) for discussions of this point). 

The plan of this section is as follows: (1) Study the equations governing the secular evolution of three coplanar Newtonian 
point masses for which Cp is small compared to the companion eccentricity and mp <^ m,, and determine the regions of 
parameter space for which the system librates or circulates. In particular, determine the fixed points of this dynamical 
system; (2) Determine how the results of (1) are modified by the inclusion of post-Newtonian terms in the potential of the 
star; (3) Introduce tidal damping and show that the system evolves towards doubly circular orbits on timescales described 
above. 



3.1 Newtonian three-body secular evolution 



Consider a copZanar point-mass three-body system with masses m«, mp and rric (for star, planet and companion respectively) 
for which rrip <ti m*, but with no constraints on rric, and let the star-planet and star-companion semimajor axes be Up and 
Oc respectivelylfl Further, let the corresponding eccentricities be and e^, with Cp <^ 1 but again with no constraint on Be, 
and let the longitudes of periastra be Wp and Wc. While no energy is exchanged between the orbits of a (nonresonant) stable 
system (and hence the semimajor axes remain constant), to first order in Bp the time-averaged rates of change of e^, ec, 
Wp and Wc are (Mardling & Lin (2002); also see iMardlind (|2007al ) for the full secular equations to all orders and including 
inclination and precession) 



dep 
lit 


15 

^ ^16 


Tip Cc 


dec 
dt 


15 / 
= ep 

16 V 


dzup 
dt 


3 

= ^np 


\ m. 


and 






dzuc 
dt 


3 

= -Uc 

4 


/rUp 



/mc\ /flp 
\m, ) \ac 



sin(t37p — ujc) 
(l-ei)5/2 ' 

sin[zup — vjc) 
(1 - e?)2 ' 



(1- 



-3/2 



) /ap\ 
\\aj 

/OpX /ep 
\ac) V ec 



COs(ci7p — ZOc) 



(l + 4eg) 
(1-e?) 



- tUc) 



(4) 
(5) 
(6) 

(7) 



where Up a nd are the orbit al frequencies (mean motions) of the planet and companion respectively. Note that lMardling fc Lin 
I 2OO2I ) and iMardlin j (|2007al ') employ a Legendre expansion in the ratio Op/ac to derive the rates of change of the elements 
and make no assumptions about the mass ratios, eccentricities or inclinations, in contrast to the literal expansion used by 
Murrav fc Dermot j jioool ) which relies on the smallness of the eccentricities and inclinations but holds for all se mimajor ax:es 
ratios . Equ ations to dT j l are given to octopole order (corresponding to the fourth power in ap/ac); see also iLee fc Peak 
( 2003h and lFord fc Rasiol I2OO0I ) . 



3.1.1 Libration versus Circulation of zup — zu^ 

In order to derive simple conditions under which the angle vjp — Wc circulates or librat e s, we assume that ec is significantly 
greater than ep ^ 1 (a situation relevant to the present study; see iNagasawa. Lin fc Idal l|2003l ) for an alternative treatment). 
Then from ((5)l, ec — 0, and we can write 



ep = — Wo ec sinr; 



(8) 



^ In fact, a Jacobian coordinate system is used to derive equations (|4]l to jTj so that the companion's orbital elements are actually 
measured with respect to the centre of mass of the star-planet system. 
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and 



n = w^-Wo 



cos?7, 



where rj = zup — zuc and Wq and Wo are constants (g for quadrupole and o for octopole) defined by 
3 ^3 



1 



acKnicJ 



and 



15 



where Cc = Vl — ei. The fixed points of (|8]) and ([9J are 

(ep,77) = (ecW^o/W^g,2fc7r) when IV, > 0, 
and 

(ep, 77) = {~ecWo/Wq, {2k + when Wq < 0. 
Here fc is any integer. Writing x = ep — Cp"^'' , where 
(5/4)(ap/ac) e^e;:^ 



ap/ac{mp/mc)ec 



and linearizing about = we have 



W 

X = -ecWo rj, 77 = ( 



a; 



^) + VP', »7 = 
which has solutions 
7j{t) = A cos{Wqt + l3), 
where 

= 7j{0f + (1 - ep(0)/e("''')' and tan/3 = (l - ep(0)/e('='')) /r;(0). 

From (|16p we conclude that the fixed point ?; = is stable. For libration around 2fc7r, k ^ 0, add 2fc7r to H17[). From 
have that 

ep(f) = 4"") [1 - A sin(M/,t + 13)] = 

SO that the maximum and minimum values of Cp when a system is librating are approximately 



max-rmn 



= (1±A) e^' 



For fibration about rj = n, also a stable fixed point, 
ri{t) =-K- A co%{Wqt + 



(9) 

(10) 

(11) 
(12) 
(13) 

(14) 

(15) 

(16) 

(17) 

(18) 
((T5]l we 

(19) 
(20) 
(21) 



with A and /3, ep(t) and e™""'"'"'" again given by JTS}, (fT9)) and (f20)) respectively. For libration around (2A: + l)7r, / 0, add 
2fc7r to H2ip . Since Wq > for libration about = and < for fibration about rj — tt, from pOf) we have that 



\ ^ l*^"- libra-tion about = 0; — ^ < t —e^ ^ for fibration about 



rj = TT. 



(22) 



Consider libration about 7; = 0. This can occur as long as there are solutions to j) = 0, that is, there exist real values for rj 
when cost; = [ep/cc) Wq/Wo- From (|19|l and 1)14^ this is (in the linear approximation) 



cos 77 : 



A sm{Wqt + P) 



(23) 



whose minimum value is 1 — yl (taking A to be positive) corresponding to the maximum value rj achieves during a libration 
cycle. But in order for ep to be positive or zero, ^ ^4 ^ 1 from (|20p . Thus we have the result that libration about 77 = 
can occur as long as \rj\ < ■k/2. Moreover, from pop the maximum value of ep associated with A = 1, and hence the value 
associated with the boundary between libration and circulation, is 
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200 




5000 



years 



Figure 1. Evolution of Cp and rj for an HD 209458 type system (parameters given in text). The full secular equations were integrated 
to produce these curves, (a) shows libratory and circulatory behaviour consistent with e^''^ = 0.063 and g^"""'*'""^ = 0.126, while (b) 
demonstrates that r] is almost constant for ep(0) not too close to g^°""'''"'y 

boundary 



The dashed curve in both plots corresponds to ep(0) 



The libration and circulation periods are consistent with the analytical value of 1.01 X 10 yr. 



boundary 



= 2e 



(=9) 



(24) 

which also holds for libration about = tt. In the latter case libration can occur as long as |?7 — 7r| < 7r/2. These results 
are consistent wi th numerically integrating the full equations (equations (|4]) to ([7]) but with the full dependence on e^; see 
Mardlin3 (2007a)), an example of which is presented in Figure[T]for a system with mp = 0.64Mj, rric = 1 Mj, m« = 1.14M0, 



Up = 0.047 AU, Oc = 0.4 AU and ec = 0.3. The values for the inner planet are appropriate for HD 209458b. Figure[TJa) plots Cp 
against for r;(0) = and 12 values of ep(0) ranging from 0.06 to 0.2. Equation ^ gives e'j^"^ = 0.063 and ej;"""'''"'^' = 0.126. 



Figure [TJb) shows that for Cp > gi'^'^ndary ^ ^ circulates with a frequency that is approximately constant (c 

boundary p„<-<-;„„. ^ _ w 



is not too close to e°™"'""'f . Putting 77 



/g gives 



77(t) = W^,t + 77(0), 
so that from ((Sj, Cp 



Wq) as long as Cp 



(25) 



-WoSin(Wgt + 77(0)) and 



ep{t) = ep(0) + (Wo/Wq) [cos{Wqt + 77(0)) - cosr?(0)] = e^,'^"-'^^ (t) 



(26) 



Therefore if Wq > 0, the maximum and minimum values of Cp correspond to = 2fc7r and 77 — {2k + l)7r respectively, where 
k is an integer, and are given by 

ep(0)±e^'=«) [lTcos77(0)], (27) 



max ,m%n 



while for Wq < 0, these are 
e^^'^^'"''" = ep(0) ± ei'"^ [1 ± cos 77(0)] 



(28) 



corresponding to 77 = {2k + 1)tt and 77 — 2kn respectively. From (|19|l and (I26p we see that both the libration and circulation 
periods are given by 2-k /Wq. For the example shown in Figure[T] 2-k /Wq = 1.01 x 10^ yr. 

The actual time variation of Ec can be obtained by forming dep/de^ from ^ and ((5|), which upon integration gives 
1/2 



eo{t) = 
where 



{I3ep{tf 



C] 



(29) 



^^1 /^/mp\ C = /3ep(0)2-v/l-e.(0)2, 

2 y ac \mc / 



(30) 

and ep{t) is given by (|19p or (|26p according to whether rj librates or circulates. 

Summarizing these results, we have that libration will occur about the fixed point (ep, 77) = {e^''\ 0) at frequency Wq > 
when 



> 



and ^ e„ < 2 el'^''' and - 7r/2 < 7? < 7r/2, 
while libration will occur about the fixed point {ep,ri) — {e^''\n) at frequency Wq < when 

and tt/2 < rj < 3tt/2. 



m, 
m. 



< 



and 



^ ep < 2 e. 



(31) 



(32) 
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CM 




Tj years 

Figure 2. Same as Figure JTJ but with Wqr included. Libration and circulation of r) is more rapid and the maximum amplitude for 
librating systems is reduced in this case. The dashed curves correspond to Cp = g^°""'''"'!' 'Yhe libration and circulation periods are 
consistent with the analytical value of 5700 yr. 



Similar results hold for fc 7^ 0. If neither of the conditions (|3H) or (|32} are met, 77 will circulate at frequency Wq with the 
maximum and minimum values of given by (f?7)l when mc/rrip > ap/ac [Wq > 0), and (pSj) when mc/m-p < ap/ac 

{Wq < 0). 



3.2 Relativistic Effects 

Contributions to the potential of the star from post-Newtonian terms affect the apsidal advance only. The orbit-averaged 
contribution to 77 from the innermost orbit is (eg. Mardling & Lin, 2002) 
3nr> ^ r, \^ 



L Gp \ C / 

where c is the speed of hght, so that 



fi = Wq + Wgr -Wo[ — ] cosn-Wq-Woi — ] cos 77, 



where 



Wq 



\acJ Vtti* 



flc V rric J 



(33) 



(34) 



(35) 



and 7 — 4:{npap/cf'(rni,/mc){ac/ap)''^ is, to first-order in ep, the ratio of Wqr to the quadrupole contribution to vjp. One 
effect of Wgr on the system is to shorten the eccentricity modulation and apsidal period (compare Figure [2] with Figure 
In this example, Wgr — Wq in contrast to Mercury's orbit around the Sun for which Wgr is_only 7% of the contributions to 



Wq from the other planets in the Solar System. The equilibrium eccentricity is now given bj[f| 



eMol\Wq\ = 



(5/4)(ap/ac 



1 — ^Jap/ac{'mp/mc)ec ^ + 'ye'i 



(36) 



Thus when apj acirapjm^e^ ^ < Ij the effect of Wgr is to decrease the equilibrium eccentricity, an example of which is 
shown in Figure [5] while for ^ apj aciviplm^eZ^ > 1 the opposite is true. Mardling fc LinI (2004) show that the former can 
be vital for the survival of short period terrestrial planets with companions. 

The analogues of (|31|l and (|32p are that libration will occur about the fixed point (ep, rj) = (ep^'' , 0) at frequency Wq > 0, 



where ep'^'^ is now given by 



when 



> \ l —£c /(I +^£c) and 



^ ep < 2 e^"«' 



and 



- 7r/2 < 77 < 7r/2, 



(37) 



while libration will occur about the fixed point {ep,r)) — (ep"^''' , tt) at frequency Wq < when 



—e^^/(l + 'yel) and ^ ep < 2 e^p"' and 7r/2 < 77 < 37r/2 



(38) 



^ Note that equation (4) in Mardling & Lin (2004) is in error by a factor of 2. 
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If neither of the conditions p7p or (|38|) are met, then when rric/mp > e^, ^ \/ ap/ac/(l+7ec), ?? will circulate at frequency in 
the positive direction with the maximum and minimum values of Cp given by (|27p . while when rric/mp < e^^ ^ apj a^j (\ + ^e\), 



77 will circulate in the negative direction at the same frequency with maximum and minimum values of ep given by (|28|l . The 
time variations of 77, ep and are given by the equations set out in the previous section, with Cp'^'^' given by (|36p. 

The form of (|36|l suggests that for a range of parameters, e^p"^^ can be large. Singular values of ep*^'' exist for some 
combinations of parameters, although when the octopole t erm in (O is included in (|9]) these are replaced by finite, although 
still relatively large values (see Nagasawa. Lin fc Ida ( 2003l l). As we will show in Section r3.3.4l this has important consequences 
for eccentricity forcing by low-mass companions. 



3.3 Tidal damping and quasi-fixed-point behaviour 

While tidal dissipation in the innermost planet directly affects the rate of change of its eccentricity, it also indirectly affects 

the rate of change of the companion planet's eccentricity (recall that the previous analysis assumed Cc was constant), as well 
as the angle between the apsidal lines, The relevant system of equations now becomes 

ep = — Wo fie sin r; — Wt Cp, (39) 

Cc = Wc ep sin rj, (40) 

ri = Wg-Wo(^^^co8rj, (41) 

where from ([2J, 

Wr^'4nJ^](^](^X-r-L (42) 



2 \ ^/p 



and from (I7|), 

W^e = i§n.(!^)(^)%-, (43) 
16 \mt./ \ac/ 

where again we have ignored the octopole term (proportional to Sp/ec) in ([7|. The damped autonomous system p9|l - (|41|l 
evolves to the formal fixed point {ep,ec,rj) = (0,0,77*), where 77* is given below and is such that |77*|(mod27r) <C or 
I77* — 7r|(mod27r) <^ for Wq > and Wq < respectively. Ifowever, there are three timescales associated with this 
evolution: (1) The system circulates until ep, ec and rj are such that libration conditions are met (see (|3ip and (|32p : the 
system is "captured" by the associated fixed point). This occurs on the circularization timescale of the inner orbit|j Tdrc', 
(2) The system librates about and evolves towards the (now quasi-)fixed point (ep, 60,77) = (ep, 6^,77*), where e^ — ec(0) and 
Bp ~ e^''^ (equation (|36|) ). on the timescale 2 Tdrc (ej! and e* will be properly defined below); (3) The system evolves smoothly 
towards the formal fixed point (cp, 60,77) — (0,0,77*) on some timescale Tc ^ Tdrc- The details for these phases are as follows, 
where for the first two we assume that ec ^ ep so that Ec — ec(0). 



3. 3. 1 Circulation 

Following Section [3. 1.11 the circulation phase can be approximated with 77 = W q. Then (|39} becomes 

ep + Wrep = -Woec sinCvygt + 77(0)) (44) 
which has solution 

ep{t) = Ce-^^* + ere*,""' cos(W,t + 77(0)) (45) 

where a = sgn \W q\ , C = ep(0) — trep'^''' cosr;(0) and we have assumed that Wt ^ \Wq\. This phase is shown in blue in 
Figure |3] and evolves on the timescale 1/Wt = Tdrc- Notice that the amplitude of variation is not afi^ected by dissipation and 
is given by e^'''^ . Also note that Qp is set artificially low in order to clearly demonstrate how this phase proceeds. 



* The contributions to ■/) from the tidal and spin bulges of the star and planets are negligible compared to the terms included here. 
^ This phase is sltipped if the initial conditions are such that the system starts librating immediately. 
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3.3.2 Libration 

The quasi-fixed point {ep,e*,Tj*) is obtained by putting = and ?) = in (|39} and (|41[) respectively. While there is only 
one true equilibrium value of Be (that is, ec = 0), we take e* to be the value of Cc when the libration amplitude of Cp drops 
below some specified small value (see next section). Eliminating ep/cc between (I39p and (|4ip then gives 

W^, \Wg\ 

where the second step comes from the fact that sgn(cos77*) = sgn(W5) and Wr/lVVq] ^ 1. Thus for Wq > 0, 77*(mod27r) ~ 
-Wr/Wq while for Wq < 0, {ri* - 7r)(mod27r) ~ -Wr/Wq. For HD 209458b with a O.lMj companion at 0.4 AU, tj* = 
10~^ — 10~^ radians, depending on the value of Cc (taking Qp = 10^). The value of Cp is given by (|36|l with Cc = e^. Linearizing 
about (e*, e*, ry*) and putting x = ep — e* and y = rj — rj* produces the system 

X = -e*Woy -Wtx, (47) 
y^{Wl/e:Wo)x, (48) 
which yields the damped librating solution 

ep{t)^e; + Ae-^'^'^^cos{Wqt + 5) (49) 
and 

rj{t) ^ri* - Ae-'^^'^''sm{Wqt + 5). (50) 

For the purposes of this analysis, the constants A and 5 can be taken to be such that the librating phase is entered at the 
first occurrence of ep — 0. Writing the time at which this occurs as t = tm,, (|49p and (|50[) become 

ep{t) = e; {l - e-^^(*-*'"')/2 cos [Wqit - tui,)]} (51) 
and 

Tj{t) - e-^^f*-*'"''/^ sin [Wq{t - tu^)] (52) 

In contrast to the circulating phase, the amplitude of oscillation of ep{t) diminishes on a timescale 2/Wt = 2rcirc as the 
system approaches the quasi- fixed point. Figure |3] shows this behaviour in red. 

For the analysis above we have assumed Ec is constant. However, as Figure |3jc) illustrates, ec slowly declines and has no 
analogue of Cp*^'' (except ec = 0). The following section studies the effect of this on the long-term evolution of the system. 



3.3.3 Long-term evolution 

Once the quasi-fixed point is reached, the companion's eccentricity evolves according to (|40p with ep = ecWo/\W q\ and 
sinr? = ~WT/\Wq\: 

WcWoWt 

■ ec 



5/2 

Wt ■ 



/25\ /mpA 
\16/ Kmc J \ac 



where 



A=(|)(^)(J)''' and F{ec) =sl{l- as:' +^ elf (54) 

with Q = ^Ja^Jocimp I rn^ (recall that Ec = \/l — ec). Note that a and 7 are effectively constant throughout the evolution 
since the semimajor axes evolve on timescales much longer than the age of the system (see Figure [5jd)). One can obtain an 
approximate timescale on which ec (and hence ep) evolves by putting Ec — sie*^) = £*, where e* is yet to be formally defined. 
Equation (|53[) then gives 

ec(t) =ec(ic)e-(*-*='/^% (55) 
where 



Tc = {F{el)/X) Tcirc 



(56) 
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years years 

Figure 3. Capture phase for HD 20 9458b with a O.lMj c ompanion at 0.4 AU with ec = 0.4. The full secular equations including 
evolution of ap and ac are integrated l lMardling fc Linll2004l ). The Q- value of the observed planet is set artificially low at Qp = 100 in 
order to illustrate the process (timescales are linearly proportional to Qp and values for e* are not very sensitive to the value of Qp; for 
this example, the relative change in e* is less than 2% over four orders of magnitude in Qp). Blue curves correspond to the circulation 
phase, red to the libration phase and green to the slow evolution or quasi-fixed point phase (see text for more discussion). In panel (c) 
notice that the relative change in is small, justifying (at least for this example) taking ~ const during this phase in the analysis. 
The diamond symbols correspond to the points (tc,e*), (tc,e*) and (tc,'?*) in panels (b), (c) and (d) respectively. Panel (d) shows the 
apsidal advance of the inner and outer orbits, with ijjp locking onto ro^. This occurs irrespective of the mass of the outer body because 
fip/fic «C 1 (see equation iTl). 



and tc corresponds to the beginning of the slow evolution phase which we choose to be the first instance of |r; — r;*] < e for 
all t > tc- Here e ^ 1 is an arbitrarily chosen parameter and is taken to be e = lO"** for the studies presented in this paper. 
This then serves as definitions of e* and Cp, that is, e* — ec{tc) and e* = ep{tc). The transition to the slow phase at time tc is 
indicated in Figures |3jb)-(d) by a diamond symbol. Note that Tc ^ Tdrc for most stable systems, except when F{ec) is small 
(see Section [3. 3. 4p . The approximate solution (|55|l is the blue dot-dashed curve in Figure |4j a). 
Keeping the full dependence on Ec in (|53p and integrating gives 



/(ec;a,7) -/(e*;a,7) = -X{t - tc)/Tc^rc, (57) 
where 

/(a;;a,7) = - ^jx^ + '^-yx^ + (a - 3^) + ^/^^ [a^ - ^7 (3x'' - llx' 23) a - i {x^ ~ 4) +'y^gix)] 

-F(7 - a + 1)^ In a; - [(7 - q)^ + l] In (^1 + (58) 

with 

g{x) = (35a;* - 185a;'' + 408x^ - 506a;^ + 563) /315. (59) 

This solution is shown as a dashed red curve in Figure SJa) and is clearly superior to the estimate (I55|l which underestimates 
the timescale on which the system evolves once apsidal locking has taken place. We define the true timescale on which 
circularization of the companion orbit takes place (the e- folding time) to be the time corresponding to e*/e, that is, from 

^true ^ [3-|-A"^(/(e*;a,7)-/(e*/e;a,7))] Tcirc, (60) 
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years years 

Figure 4. Continuation of numerical solution in Figure |3] showing the long-term evolution of (a) Cc and (b) ep (black curves, sitting 
underneath rod curves), again with Qp = 100 for illustration (note the y-axis scales are different to those in Figure [S]!. This corresponds 
to 125 Gyr of evolution when Qp = 10^ (compare this to Tdrc = 0.045 Gyr). Also plotted in (a) is approximate solution I I55I I (blue 
dot-dash curve) and the analytic solution l|58|l (dashed red curve). The approximate solution underestimates the timescale on which ec 
evolves by a factor of 1.8. Also plotted in (b) is the analytic expression for ep{t) (red curve) using II36I I with Cc = ec(t) from II58I I. The 
analytic solutions are essentially indistinguishable from the numerical solutions. 



where we have put tc — STdrc (1 for the circulation phase + 2 for the libration phase). As ec slowly decreases, Cp maintains a 
value given by (|36|l . Any departure from this quasi-equilibrium value is quickly damped out on the circularization timescale 
Tcirc- Thus we can use (|57|l and (|36p to calculate analytically ec{t) and ep{t) during the slow phase, and hence to determine 
the current value of Bp using the estimated age of the system. The analytic solution for ep{t) is shown in Figure |4jb) (red 
curve) and is indistinguishable from the numerical solution. 



3.3.4 The increasing ep regime: very low-mass companions 

The solution for the example given in Figure |4] is essentially a modified decaying exponential. However, there is another 
regime of the function f{x; a, 7) which behaves very differently and is associated with very low-mass companions which, for 
a large part of the life of such systems, force Cp to increase to some maximum value before decaying to zero at the end of 
the tidally evolving phase. An example of this is shown in Figure [5] for the HD 209458 system with a companion of mass 
rUc = 0.0125Mj — 4:.2M^ at = 0.3 AU and an initial companion eccentricity of 0.7, a high value chosen to clearly illustrate 
the theory. The system enters the slow phase with e* = 0.705, e* = 0.013 and rj* ~ — 437r. The average reflex velocity of the 
star due to the companion is 0.64 ms-^ See the fi gure caption for a detailed description of the evolution of the system. The 
red dashed curves again represent the analytic solutions (|58|l and (|36p . clearly demonstrating their predictive power. In fact, 
the analytic solution for Cp passes through a singularity as the denominator in (|36|l passes through zero. The full numerical 
solution avoids this, but still manages to increase to a signiflcant level on a reasonably long timescale. For Qp — 10^, passage 
through the "singularity" occurs at 2.2 Gyr, well before the system age of 5.5 Gyr, while for Qp — 2x 10^, passage through this 
point occurs "now" . Note that this system is stable when the full three-body (non- secular) equations are integrated including 



post-Newtonian relativistic terms, tidal and spin bulges, as well as tidal damping ([Mardling fc Linll2002[ ). 

Although systems which pass through the singularity suffer a rapid decline in Cp once they have done so, they tend to 
maintain the value of e* that they had at the beginning of the slow phase before this happens (see also Figure |6]). This is in 
contrast to systems which don't pass through the singularity; these simply decay approximately exponentially. To determine 
the parts of parameter space for which Cp necessarily passes through the singularity and therefore maintains before it does 
so, and in particular which systems pass through the singularity after the current age of the system, we start by studying the 
behaviour of the denominator of (|36|l . Defining 

A(ec;a,7) = l-aec~'+7£^ (61) 



where (again), a = ap/ac{mp/mc), Figure [7| plots this function against rric for specific values of Oc and for Cc = 0.8 and 
Be = 0. There are three regimes of interest. As the system evolves from e* down to Cc — 0, either (a) A (and hence Wq) 
remains negative so that cost; ~ — 1 always (blue arrow); (b) A passes through zero so that rj goes through a short period 
of circulation as cost; passes from -1 to 1 (yellow background, red arrows) or (c) A is positive always and cos?; ~ 1 (green 
arrows). In the case shown here (e* = 0.8), systems cross A = for a finite range of values of rric when Oc 0.29; when 
0.29 ^ Oc < 0.52, all systems with rric less than some critical value cross A = 0; when Oc >, 0.52 no systems cross A = 0. 
Figure [Tjb), for example, shows that all systems with nic <, 0.25i\fj and Oc — 0.3 (and e* = 0.8) necessarily evolve through 
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Figure 5. A terrestrial- mass planet can sustain the observed eccentricity of HD 209458 - but not for long enough in this example where 
rUc = 4.2M0, at least for Qp 2 X 10^. Recall that the actual timescale is proportional to Qp so that 6 Myr with Qp = 100 corresponds 
to 6 Gyr with Qp = 10^. The black curves show the numerical solutions while the red dashed curves are given by the analytical solutions 
l|36|l for Bp and Il58|l for Be- (a): The evolution of Cp following capture. After effectively maintaining the value e* ~ 0.013 for some 2 
Gyr (for Qp = 10^), ep increases to a maximum of around 0.04 before rapidly declining to around 0.008 and then slowly declining to 
zero. The increase in ep corresponds to a singularity in the analytic solution which fairly accurately follows the full solution. The fact 
that the analytical solution is reliable confirms this dependence. Thus, for Qp = 10^, the maximum value of ep occurs around 2.7 Gyr. 
Recall also that the value of Bp at the beginning of the slow phase is quite insensitive to ep(0), 77(0) and Qp. (b): The evolution of 
following capture. Again the analytic solution is a good approximation to the full solution. The blue dot-dashed curve is the approximate 
analytic solution given by I I55I I. Thus when a system must pass through the singularity, the timescale given by I I56I I overestimates the 
circularization timescale in contrast to the previous example, (c): The evolution of r; = zup — zu^. During the pre-capture phase, rj 
circulates backwards before being captured at = — 437r (recall that when Wq < 0, cosry ~ — !)• As the system passes through the 
singularity, that is, as Wq passes through zero to become positive (see equation ||35}), rj circulates back up to — 407r so that cos 77 ~ 1. 
The brief circulatory period is accompanied by oscillations in Bp and Be which can be seen in (a) and (b). (d): The evolution of Op. The 
full (secular) numerical solution includes the evolution of both semimajor axes; here we see that over the entire evolution of this system. 
Op changes by less than 4%. 



A = 0. The ranges above vary with different values of e*; curves corresponding to e* < 0.8 he between the two shown. What 
remains to be determined, however, is which of these systems maintain a non-zero value of Cp for at least as long as the current 
age of the system before they do so. Onc e this is done, the parameter ranges should be further restricted by determining 
which systems are stable (|Mardling||2007br ). 

Figure |8] plots ep(now), ep{tc) and ec(now) against companion mass for the HD 209458 system for various values of Uc and 
e*. Here "now" refers to t = 5.5 Gyr, and the quantities are calculated using (|36|l and (I58|l . The green-filled circle corresponds 
to the system shown in Figure \Ba.t t = 5.5 Gyr for Qp = 10^ {t = 5.5 Myr for Qp = 10^), while the red-filled circles 
correspond to Figure [6] at the same time. Note the extremely low companion masses in Figures |6][a) and (b): rric = I.9M0 
and rric = O.4M0 respectively. The maximum values of ep achieved are (a) 0.022 at 5.8 Gyr and (b) 0.012 at 6.0 Gyr. The 
corresponding values of Cc are 0.27 and 0.50. At 5.5 Gyr, the values for ep are (a) 0.010 and (b) 0.004, with corresponding 
values for ec of 0.41 and 0.60. Thus scenario (a) 
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HD 209458: 
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Figure 6. Similar to Figure[5]but with lower values for rric. Passage through A = occurs later in the life of the system (compare panel 
(a) to Figure [Sfa)), while e* is smaller. For (a), the maximum value of Cp = 0.022 is reached at 5.8 Gyr for Qp = 10^ (ie., 5.8 Myr for 
Qp = 100), while at 5.5 Gyr, Cp = 0.010. For (b), Cp = 0.012 at 5.96 Gyr, while ep = 0.004 at 5.5 Gyr. 
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Figure 7. Behaviour of the denominator of l|36| l. A, used to determine which configurations have e^''"' necessarily passing through a 
"singularity" during their evolution. See text for discussion. Note tha t some of th e systems represented here are not stable when the full 
three-body equations are integrated with appropriate damping (iMardling &: Linir2002i1 . 



while scenario (b) requires a Q-value around 10% smaller in order that the maximum value of Cp occurs about now (recall 
that this also depends on the age estimate for the star). Figure |8] demonstrates that high initial values of e* are required in 
order that extremely low-mass companions be responsible for maintaining the eccentricity of inflated planets (see Section [6] 
for a discussion on this point). 

Figure |8] together with Figures [5] and [6] and equation (|60|) , suggest that for given values of Oc and e* , the time a system 
takes to achieve its maximum value of ep (given it passes through A = 0) increases as the companion mass decreases. Figure[9] 
plots this time, Tcvoi, as a function of rric for various values of Oc and (solid curves). Also plotted is r^'^'"^ (equation (|60|l ) 
which is the e-folding time for systems which don't go through A = but simply decay approximately exponentially (dashed 
curves). The small values of Tcvoi corresponding to, for example, systems for which e* — 0.5 and flc = 0.34 with rUc <, 0.002Mj 
(the dark blue curve which appears coincident with the a;-axis), compared with the large values of Tevoi for rric ^ 0.003A/j, 
simply indicate that passage through A = occurs early in the life of the system, while decay following passage is very slow. 
Figure [10] illustrates this with the above example. The blue diamond indicates the nominal beginning of the slow phase (as 
defined in the discussion following equation (I56p ). There is little difference between Figures [TOT a) and (b) except that lior bl 
does not meet the slow phase criterion until after the system passes through A = 0. 

It is of interest to ask whether or not it is possible for a system to skip the librating phase and evolve directly to the true 
fixed point (ep,ec) — (0,0). Referring to equations (|37p and (|38|l we see that for Cp < 2ep'^''' a system will librate as long as 
— 7r/2 < rj < n/2 or n/2 < rj < 3tt/2. Since a circulating system must pass through all values of 77, once ep is permanently 
below 2ep'^''' it has no choice but to librate until a local quasi-fixed point is reached. 
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HD 209458: age = 5.5 Gyr, Q^, = 10^ 
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Figure 8. The values of Cp and ec (black solid and red-dash curves respectively) at t = 5.5 Gyr, and the value of ep at the beginning 
of the slow phase (ep(tc); blue dot-dash curves), all with Qp = 10^. From Figures |5] and [6] it is evident that systems which pass through 
A = in their lifetime approximately maintain their value of ep at 4 = before they do so. Thus, for example, in (a) all systems with 
rric < 0.004Mj will have the same finite (but small) eccentricities they had when they entered the slow phase. For a small range of 
masses (0.004 < m^/Mj < 0.006) ep(now) > e* (the black curve sits above the blue curve), while for greater masses ep(now) ^ e*. Note 
that formally these curves have no upper bound at the point corresponding to A = 0, however the finite resolution of the plots suggest 
erroneously that they do. In reality there does exist an upper bound when the secular equations are directly integrated as is evident 
in Figures [5] and |6] It is clear from this figure that high initial values of e* are required to have very low-mass companions responsible 
for maintaining the eccentricity of inflated planets, although lower values of e* are required for lower values of ac (see Figure 1131 and 
Section |6] for a discussion on this point). 



4 CONSTRAINTS ON ORBITAL PARAMETERS OF COMPANIONS 

We are now in the position to constrain the orbital parameters and masses of possible companion planets to HD 209458b, 
HAT-Plb and WASP-lb. Given estimated ages, our aim is to find initial conditions and companion masses for which these 
planets have eccentricities in the ranges observed (WASP-lb is currently unconstrained), and in the case of HD 209458, 
produces a stellar reflex velocity l ess tha n 10ms~^ (an upper limit which gives companion masses less than the upper bound 



of 0.3Mj cited in lLauehlin et al.l (|2005aD). 

Acknowledging again that these systems are unlikely to be exactly coplanar (but see Figure [15] which shows that moderate 
inclination has little efltect), we assume zero mutual inclination and vary the parameters ac, rric and ec(0). Our results 
correspond to the case ep(0) = 0.1 and rj — 170°, but are not very sensitive to using other values of these parameters (ie, they 
evolve to the quasi-fixed point {ep,ec,rj*) initially with a variation in e* of at most around 10%). Note also that while our 
aim is to determine whether or not current-day solutions exist (especially for HD 209458 which is the most constrained), we 
can also say something about the orbital elements at earlier times and hence (cautiously) constrain the initial conditions. 
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Figure 9. The time to reacli tiie maximum value of ep, Tg^oj, for systems which pass through A = (solid curves), and the evolution 
timescale (defined as the time corresponding to ec = e*/e) for systems which do not (dashed curves). Various values of ac are indicated 
on the plots; the estimated age of HD 209458 is indicated by the horizontal dashed line. The lower the companion mass, the longer the 
system takes before it passes through A = (given that it does so at all). This phenomenon favours systems with high values of e*. 
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Figure 10. Demonstration that for some systems, passage through A = occurs very early on in the evolution, and circularization is 
approximately exponential after that. The blue diamond symbols again indicate the formal beginning of the slow phase. 



4.1 HD 209458 

Figure [TTI presents results for 7?(0) = 170°, 6^(0) = 0.1, 0.3 < ec(0) < 0.8, 0.02 ^ m^/Mj ^ 0.2 and 0.1 < aJAU < 0.8. 
The full secular equations 

used in (|58p to calculate the current value of ec ("fie (now)" ) and hence of ("ep(now)") via p6|) for a given actual Q-value 
{Qp ~ 10^). Recall from the discussion in Section [33] that timescales scale linearly with Qp and that the relative change m 
is less than 2% over four orders of magnitude in Qp. Each point in the rUc — flc plane is coloured according to epijnow) (top 
six panels) and Cc (now) (bottom six panels) . Also shown are boundaries for stellar reflex velocities at periastron — 10ms~^ 
and = 5ms'^. Defining a "solution" to be a conflguration for which ep(now) > 0.010, solutions exist below — lOms"^ 
for ec(0) > 0.3, ec(now) > 0.2, Oc > 0.35 AU and > 0.05 Mj, and below v^, = 5ms"^ for 6^(0) ^ 0.5, ec(now) > 0.3, 
ac >, 0.4 AU and rUc >, 0.05 Mj. However, notice the orange patches around ac=0.2 AU corresponding to nic < 0.03 Mj: 
enlargements of these are shown in Figure \TI\ The slow phase of the evolution of these systems actually involves an increase 
in Cp as discussed in Section 13.3.41 so that companion masses as low as IMg are capable of exciting signiflcant eccentricity 
in the observed planet. Some systems reach eccentricities around 0.02 (the aqua points in the top three plots of Figure [T2|l . 
This effect is more pronounced when Qp is higher as is demonstrated in Figure [13] for the case Q = 3 x 10^. 



Mardling 2007a ) were integrated with Qp 



100 until e"! was reached and the value for this was 
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Figure 11. "Current" values of the eccentricity of HD 209458b, ep(now), (top three panels) and of its hypothetical companion, ec(now), 
(bottom three panels) for a range of hypothetical companion masses and semimajor axes. An estimated age for the system of 5.5 Gyr 
is assumed as well as Qp = 10®. Initial values for are indicated at the top of each panel, while ep(0) = 0.1 and r;(0) = 170° for all 
experiments. The minimum companion mass is rric = 0.02Mj = 6.67A^g3. Note the circle on the panels corresponding to ec(0) = 0.7: 
although the companion orbit starts with a high eccentricity, its current eccentricity is relatively low (around 0.35). Also shown are 
curves of constant stellar reflex velocity, rij, due to the companion at periastron. See text for discussion. 



4.2 HAT-PI 

HAT-PI is a member of a binary system for which the projected separation of the stars is 1550 AU ( Bakos et al.ll200^ . At this 



distance, the companion is incapable of inducing any eccentricity in the planet's orbit because the tidal damping timescale is 
far shorter than any timescale on which eccentricity could be induced. 

As far as the author is aware, at the time of writing no upper bound has been put on the refle x velocity of the st ar due 



to a companion planet. The preliminary value of Bp = 0.09 ± 0.02 awaits reflnement, however, as iBakos et al.l lj200d ) note, 
it produces a similar tidal dissipation rate in the planet to HD 209458b, consistent with its similar radius and its increased 
distance from the star (see Table[2}. We have performed similar experiments for this system as was done for HD 209458, except 
for a wider range of companion masses and semimajor axes. The results of these are presented in Figure [14] for Qp — 10^ . 
These results suggest that if a companion planet is responsible for the inflated radius of HAT-Plb, it should be detectable 
in the radial velocity signal, except if a very low-mass companion is responsible, solutions for which are similar to those for 
HD 209458 (not visible in Figure [14]). If a companion planet is detected and a reliable measurement of the eccentricity of 
HAT-Plb is made (as well as reliable estimates for the companion's mass, semimajor axis and eccentricity), it may be possible 
to use the theory developed here to deduce a relationship between the Q-value of HAT-Plb and e*, and moreover, to put a 
lower bound on the value of Qp (recall that e* is the companion's eccentricity when the system enters the slow phase and 
hence is an an estimate of the initial eccentricity). The procedure is as follows. 



(i) Check how accurately equation (|36|l giving the relationship between and ap, ac, rup, rric and Cc holds. Any substantial 
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Figure 12. Detail of the lower left-hand corners of two of the panels in Figure [TTI ieriff) = 0.70) together with those for ec(0) = 0.60 
and ec(0) = 0.80. These indicate that "solutions" exist even for companion masses as low as an Earth mass, and that the number of 
solutions increases with increasing ec(0). Note in particular the blue dots in the top panels indicating relatively high eccentricities: these 
correspond to the maxima see in Figures Isl and l6l The solution circled corresponds to rric ~ I.OMq with ep(now) = 0.03, ec(now) ~ 0.5 
and ac = 0.3 4 AU. Note also that not all systems represented here would be stable had a direct integration code been used to do the 
integrations l|Mardling fc Linll2004h . The minimum companion mass is rric = O.OOlMj = O.33M0. Also shown are curves of constant 
stellar reflex velocity, due to the companion at periastron. 



departure will be due to either (a) a nearby third planet; or (b) a significantly non-zero mutual inclination; or (c) a Q-value 
which is large enough for the system not to have yet entered the slow phase. In case (a), a third planet will modulate the 
eccentricity of its inner siblings, thereby preventing the equilibriating process from occurring (unless 3r*''"'^ is less than the 
age of the system, in which case the the process can work on all three planets!) In case (b), the effect of mutual inclination 
is shown in Figure 1151 The value of e* is slightly reduced, and the eccentricity is no longer non-oscillatory during the slow 
phase, although the oscillation amplitude is small for moderate mutual inclinations. The mutual inclination itself oscillates 
with an amplitude of around 20°. Inclinations need to be greater than around 40° to signficantly affect the evolution. In 
case (c), since tc — STdrc — 0.33 Gyr for Qp = 10^, Qp would have to be greater than 10^ to still be substantially far from 
pseudo-equilibrium (given an age of 3.6 Gyr). 

(ii) Given that (I36|l is satisfied, (|57[) can be used with t — 3.6 Gyr and tc — 3 Tdrc to give a relationship between Qp and 
e* and in particular, to put a lower bound on Qp. Examples of this are plotted in Figure [16] for rric — IMj and for several 
values of ac and ec(now). For example, if a companion is found at Oc — Q.4AU with Be = 0.5 (corresponding to Cp = 0.11), 
one can conclude that Qp > 1.0 x 10^. 
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HD 209458: age = 5.5Gyr, 2 = 3x10 
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Figure 13. Similar to Figure [121 but witli Qp = 3 X 10'' sliowing tiie dramatic effect increasing Qp has on ep(now). Note the general 
reductio n in the values of Be from edO) to ee(no w). Solutions exist for ec(0) ^ 0.4, an initial eccentricity not inconsistent with those 
found bv lMandell. Raymond. &: SigurdssonI ||2007| ): higher values of ec(0) yield more solutions. Also shown are curves of constant stellar 
reflex velocity, wj, due to the companion at periastron. 



4.3 WASP-1 



While iBurrows et al.l (|2007|) find that around four times as much energy is required to inflate WASP- lb to its observed size 
than for HD 209458b (they used a value for its mass of 0.87i\fj; its revised estimate is 0.79Mj compared to 0.64A/j for HD 
209458), its closer proximity to the parent star results in a similar value for Cp for both systems (see Table [T]). While there is 
currently no published estimate for the eccentricity of WASP- lb, numerical experiments performed for this system with an 
estimated age of 2 Gyr show a range of solutions very similar to those for HAT-PI (Figure [14)) but with the legend of values 
for ep(now) replaced by one with a maximum of 0.05 instead of 0.1. Very low-mass solutions also exist for this system, albeit 
for a smaller range of parameters (or higher Q- values) than for HD 209458b and HAT-Plb. Note that Table [2] suggests Cp is 
likely to be similar to that of HD 209458 if a companion is responsible for its inflated radius. 

In contrast to HAT-PI, it almost certainly will be possible to put a lower bound on the value of Qp if a companion 
planet is detected and a reliable measurement of the eccentricity of WASP- lb is made (as well as reliable estimates for the 
companion's mass, semimajor axis and eccentricity). This is because Tdrc is only 14 Myr for Qp — 10^, so that Qp would 
have to be greater than 2.4 x 10^ for the system to be substantially far from the equilibrium phase. Figure [17] is similar to 
Figure [161 however note the increase in Qp for the case ec(now)=0.5, flc = 0.3AU. This is consistent with the fact that Tdra is 
relatively short, and since r*'^"'^ (equation ((60|) is proportional to Tdrc, close companions will also circularize rather quickly 
unless Qp is large. 
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Figure 14. "Current" values of the eccentricity of HAT-Plb, ep(now), (top three panels) and of its hypothetical companion, ec(now), 
(bottom three panels) for a range of hypothetical companion masses and semimajor axes. An estimated age for the system of 3.6 Gyr 
is assumed as well as Qp = 10^. Initial values for are indicated at the top of each panel, while ep(0) = 0.1 and r)(0) = 170° for all 
experiments. Initial values for ec are indicated at the top of each panel, while ep(0) = 0.1 and 77(0) = 170° for all experiments. Note the 
scale used for ep(now) is different to that used for HD 209458. The tentative observed eccentricity is 0.09 ± 0.02. Also shown are curves 
of constant stellar reflex velocity, due to the companion at periastron. See text for discussion. 
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Figure 15. The effect of inclination. When the orbits are mutually inclined the value of e* is reduced slightly, and the eccentricity is 
no longer non-oscillatory during the slow phase, although the oscillation amplitude is small for moderate mutual inclinations ( < 35°). 
Kozai cycles are present here for v^j = 40°. This example is for a HAT-Pl-type system with ac(0) = 0.4, ec(0) = 0.5, rric = O.bMj and 
Qp again set artificially low at 100. 
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Figure 16. Using observed data to put a lower bound on Qp for HAT-Plb. Example shown here is for rric = i-Mj. See text for discussion. 
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Figure 17. Using observed data to put a lower bound on Qp for WASP-lb. Example shown here is for rric = IMj. See text for discussion. 



5 CIRCULARIZING THE ORBITS OF PLANETS IN THE HABITABLE ZONE 

Motiv ated by the work of Raymond. Mandell fc Sigurdsson ( 2006h , Mandell. Raymond, fc Sigurdsson ( 2007 ) and Fogg fc Nelson 
1 2006l l who find Earth- mass planets forming externally to a migrating gas giant after it scatters planetesimals outward past 
its own orbit, we use the analysis presented here to identify two-planet configurations for which the outer planet resides in 
the habitable zone of the system, and for which the circularization timescale of that planet would be longer than the age of 
the system without the assistance of a short-period companion. Using equation (|57|) with Cc — 0.02 (it is singular at Be — 0), 
one can determine which values of rric/mp and a^jap produce t < tage, where tage is the age of the star. Figure [181 shows an 
example of this for the case Up — 0.04, mp = 0.6Mj, = 0.3 and m* = 0.6Mq. These parameters were chosen because the set 
of companion masses and semimajor axes for which the companion has already circularized contains a subset of Earth-mass 
planets which lie inside the habitable zone of the star, defined to be such that the semimajor axis lies in the range 

0.8 (L./L©)'''' < ac/AU < 1.5 (L./Lq)'/' , (62) 

where L* is the luminosity of the star and we have followed the definition given in IScalo et all l|2007 l) which gives an empirical 
relationship between the luminosity and the mass of a star. The habitable zone also contains a substantial fraction of planets 
with masses up to IMj; these systems are also of interest because they can harbour moons which could be habitable. Figure [T^ 
shows the effect of decreasing and increasing the stellar mass. Increasing e* pushes the minimum of the circularized region 
upwards; increasing the mass of the short-period planet pushes this region to the right as does increasing ap. 

Currently, there are no known isolated short-period pairs of planets to which this theory can be applied. However, 
low-mass stars are now being targeted more frequently so it seems likely such systems will be discovered in the near future. 
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Figure 18. Circularized planets in the habitable zone. If the condition t < t^ge is met, a coloured dot is plotted according to the scheme: 
red for tags = 6 Gyr, green for for tage = 3 Gyr, and blue for for tage = 1 Gyr. A Q-value of 10^ and a radius of 1.3-Rj are assumed for 
the short-period planet. The hatched region corresponds to the habitable zone of the star. For these systems, ap = 0.04, nip = O.OMj, 
e* = 0.3 and m* = O.GMq. Panel (b) shows detail of panel (a) for low-mass planets {M^ is the mass of the Earth). 
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Figure 19. Same parameters as in Figure [TSl but with (a) m* = O.^Mq and (b) m,» = O.TA/q. 



6 SUMMARY AND DISCUSSION 
6.1 Summary 

The ideas and results presented in this paper are summarized as follows: 

(i) If a planet is inflated because it is tidally heated, its radius is related to the orbital eccentricity by equation ([3}. The 
biggest uncertainties are the Q-value of the planet, Qp, and the heating needed to maintain the planet's size, Cp. The latter 
depends on the structural model used, and in particular, where inside the planet the energy is dissipated. 

(ii) Conditions for libration and circulation of r; = vjp — zoc including relativistic effects are given by (I37|l and (|38p for the 
case that ep <^ Ec and nip ^ m, (there are no conditions on rric), with the time dependence of Cp, Cc and 77 given by (|19p . 
([261). P9ll. ini), (HI]), and ((25|, with e^'^"' given by (|36ll. 

(iii) The system evolves through three stages when tidal dissipation is present (Figure [3]): (1) circulation of the angle 
between the lines of apsides of the two planets accompanied by the slow oscillation of the eccentricities at constant amplitude, 
together with decline of the mean value of the inner planet's eccentricity until it reaches a (quasi)-fixed value (equation (145^ ): 
this occurs on the circularization timescale; (2) libration of the angle between the lines of apsides accompanied by the slow 
oscillation of the eccentricities with reducing amplitude but maintenance of the mean value of the inner eccentricity (equations 
(|49l) and (|50p with the mean value ep'^''' given by (|36j): this occurs on twice the circularization timescale; and (3) either (a): 
a slow non-oscillatory decline in both eccentricities to zero or (b): approximate maintenance of the inner eccentricity at the 
constant value ep = e* for some (often extremely long) period of time followed by its increase to some maximum value followed 
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by an often fast decline of both eccentricities to zero. Both are described by (|36p together with (|57[). Whether or not the 
final phase occurs via route (a) or (b), it proceeds on a timescale sometimes several orders of magnitude longer than the tidal 
circularization timescale of the inner planet, Tdrc- 

(iv) As a consequence of the existence of route (b) above, very low-mass companions are sometimes capable of sustaining 
the eccentricity needed to inflate short-period gas giants (eg. Figures [5] and |6} . For this to occur it is at least necessary that 
the companion eccentricities are sufficiently high; the lower the companion mass the higher the initial eccentricity needed. 
The timescale on which the eccentricity is sustained increases with decreasing companion mass, however, the value of the 
sustained eccentricity decreases with dec reasing mass. For a gi ven s ystem age, more solutions exist fo r high er values of Qp. 
This scenario is supported by the work of Fogg fc Nelson ( 20061 ) and Mandell. Ravmond. fc Sigurdsson (2007) (see discussion 
below) . 

(v) As well as the very low-mass solutions produced by route (b), route (a) produces solutions for which the eccentricity 
and hence radius of HD 209458b is sustained by a companion which produces a stellar reflex velocity below the stellar jitter 
estimate (Figure [11}. More solutions exist for higher Q-values of HD 209458b. 

(vi) Companion planet solutions are presented for HAT-PI and WASP-1 for a wide range of companion parameters. 

(vii) If companions to HAT-Plb and/or WASP-lb (or other systems like these) are discovered and reliable estimates of the 
system parameters can be made (and there is no sign of any other companions), it should be possible to put lower bounds on 
the Q- value of the short-period planet in the system (Figures 1161 and fT7|) . Moderately inclined orbits have little effect on the 
long-term evolution. 

(viii) Long after the protoplanetary disk has disappeared, it is possible to circularize the orbits of some planets whose tidal 
circularization timescales are longer than the age of the system if they have short-period companions whose own values of 
Tcirc are short. The circularization timescale of a companion is given by equation (|60p . and depends on the structure of the 
short-period planet. 

(ix) In the light of work bv lMandell. Ravmond. fc SigurdssonI 1 20071 ) and [Fogg, fc NelsonI |2006h . systems with hot Jupiters 
orbiting stars with masses around 0.5 — O.7M0 are good candidates for finding circular Earth-mass planets in the habitable 
zone, as long as there is no third planet nearby. 



6.2 Discussion 



How much the early evolution of short-period two-planet systems proceeds according to t he simple picture described here 
depe nds on how planets come to be in their current positions. If a gaseous disk pilots them in (ILin. Bodenheimer fc Richardson 
1996h it will influence this phase at some level, especially in a case like WASP-lb for which T drc is as low as 14 M yr (a time 
comparable to the disk lifetime). If, instead, a planet is scattered in to its current position ( Rasio fc FordI 1996h it is likely 



to arrive with significant eccentricity which must subsequently be damped. Recall that the analysis here assumes Bp ^ Cc, 
although relaxing this condition makes little difference to the overall picture. 

In addition, conditions at the time of arrival will be very different to those in a mature system. For example, a planet's 
radius will be significantly larger, thereby increasing the strength of the tidal interaction with the star (see Fig. 1 in Bo- 
denheimer, Lin & Mardling 2001 for the time-dependence of the planetary radius for several models). The stellar spin rate 
will be significantly higher, and in some cases be high enough to force an increase in the planet's eccentricity. The lat- 
ter may account for the range of ec centricities observed in the "borderline planets" with periods in the range 7-21 days 
( Dobbs-Dixon. Lin fc Mardling! [iooi ). 

It s eems likely that Nature uses a hybrid of migration and scattering, thereby producing a rich variety of possible configu- 
ratio ns (iMooiheadfc Adamsllioosi : IPogg fc Nelsonlbood : iMandell. Ravmond. fc Sigurdssonll2007l : iNagasawa. Lin fc ThommesI 



2OO5I : IZhou et al 



Of particular interest here is the situation where planetesimals are shepherded in by a migrating 



gas giant, these going on either to accum ulate to form a hot Earth internal to the gas giant's orbit, or to be scattered out 
before accumulating in an external orbit ( Mandell. Ravmond. fc Sigurdsson 2007 : Fogg fc Nelson 20061 ). In the former case 
the hot Earth is often subsequently scattered out, so that both scenarios produce eccentric low-mass objects with periods a 
few times that of the massive planet. The final eccentricities of these Earth-mass planets tend to be fairly low due to gas 
drag and dynamical friction with nearby planetesimals, however the simulations of the above authors sampled only a small 
range of initial condi tions and it may be that substantially hi gher eccentricities can result from this process. The maximum 
eccentricity found bv lMandell. Ravmond. fc SigurdssonI 1120071 ) was 0.46 for a planet with mass 0.32M® at Oc = 1.74 AU out- 
side a Jupiter-mass giant at Op = 0.14 AU (simulation JD-3 in their Table 4). Higher eccentricities may result when low-mass 
planets are scattered in via encounters with bodies further out in the disk. Note that the semimajor axis ratio in the above 
example was 12.4, about twice that of the Earth-mass "solutions" for HD 209458 (Figure [13]) , while others were substantially 
lower. 

While extremely low-mass solutions are relatively i mprobable, Fig ure ll3l indica tes that there are many low-mass solution s 
for HD 209458 which are consistent with the findings of lFogg fc Nelso n (2006) and Mandell Ravmond. fc Sigurdsson (12007 ), 
especially for such Q-values (3 x 10^). While the findings of Bodenheimer. Laughlin fc Lin ( 20031 ) and Burrows et al. (2007) 
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suggest that only low Q-values are consistent with the observed eccentricity (if indeed it is non-zero) , it should be remembered 
that these authors distributed the tidal energy homogeneously throughout the planet's envelope and implicitly assumed that 
the energy was dissipated locally. If, in re ality^ the energy is dissipated near the surface (but not so close to the surface that 
there is essentially no structural change i Wul 20051 ) ). then less tidal energy may be needed to produce the observed radius 
for a given Q- value. This allows for higher Q- values to be associated with a given value of ep so that solutions such as those 
shown in Figure [131 mav. in fact, be viable. If so, inflated planets are good targets in the search for terrestrial planets. 

Of course, while very low mass solutions are of extreme interest (and completely unexpected!), there are still many 
higher companion mass solutions for HD 209458b, especially for Q-values above 10^. While WASP-lb and HAT-Plb remain 
unconstrained (except for the latter's eccentricity), there exist almost unlimited potential solutions all of which offer the 
opportunity to put lower bounds on their Q-values if companion planets are discovered. Moreover, if companions are identified 
and their orbital parameters accurately measured, the analysis presented here can be used to constrain their history, in 
particular, their "initial" eccentricity. 

Conclusions such as these are possible as long as there is no third planet nearby, and as long as one is confident that the 
circularization timescale (or more precisely, three times the circularization timescale) is shorter than the age of the system. 
Another possible departure from the assumptions made in this analysis is non-coplanarity. However, unless the relative 
inclination of the two planets is substantial, say, greater, than around 40°, non-coplanarity makes very little difference to the 
predicted outcome, as has been verified by integration of the full secular equations (Figure [Tsjl . 

Remarkably, at this stage there are no known two-planet systems for which the innermost planet has a circularization 
timescale less than a Gyr, and thus it is currently not possible to compare our two-planet model with a real system. In the 
meantime, equations (|36p and (|58p provide clear guidance on how such systems depend on the system parameters. 
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